Localization of Light: Dual Symmetry between Absorption and Amplification 
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We study the propagation of radiation through a disordered 
waveguide with a complex dielectric constant e, and show that 
dual systems, which differ only in the sign of the imaginary 
part of e, have the same localization length. Paradoxically, 
absorption and stimulated emission of radiation suppress the 
transmittance of the waveguide in the same way. 
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Localization of electromagnetic waves in a random 
medium has attracted much interest since the origi- 
nal proposals of John (2) and Anderson jjj. An essential 
difference with localization of electrons is the absence 
of a conservation law for photons. Light is absorbed or 
amplified — while retaining the phase coherence — if the 
dielectric constant has a non-zero imaginary part. The 
intensity of the radiation which has propagated without 
reflection over a distance L is then multiplied by a factor 
e aL , with a negative (positive) for absorption (amplifi- 
cation) . The interplay of absorption and localization has 
been studied extensively For the one-dimensional 

problem of a disordered single-mode waveguide (length 
L, mean free path I), the result for the transmittance T 
(being the ratio of transmitted and incident flux) is 



(InT) = (a- l~ x )L, 



(1) 



where (• • •) denotes an average over disorder. Eq. (|lj) was 
derived for a < 0, corresponding to absorption. 

In this paper we address the question: What happens 
for amplification? Since T = e aL in the absence of reflec- 
tion for both positive and negative a, one might surmise 
that Eq. ([j]) holds both for absorption and amplification. 
This is correct for short waveguides. However, we find 
that the asymptotic result for L — > 00 is 



(lnT) = (-H-r 1 )L + 0(l). 



(2) 



We will show that exponential decay of the transmittance 
in the case of amplification, (InT) ~ —L/£, is in fact 
implied by its exponential decay in the case of absorption, 
with a duality relation between decay lengths: 



£(a)=£(-a). 



(3) 



This duality relation extends beyond the strictly one- 
dimensional case, the only essential ingredient being an 



exponentially decaying transmittance in an absorbing 
system. Contrary to intuition, amplification suppresses 
the transmittance in the large-L limit just as much as 
absorption does. 

Experimentally, a random amplifying medium can be 
realized in a turbid laser dye or a powdered laser crystal 
. Stimulated emission of radiation leads to a dielec- 
tric constant with a negative imaginary part, correspond- 
ing to a > 0. We do not present a complete theory for 
such a "random laser" , because we ignore spontaneous 
emission. (This would correspond to a source term in 
the wave equation fl2|| , which we do not include.) Still, 
because of the different time scales for stimulated and 
spontaneous emissions, we believe that a time-resolved 
experiment in a waveguide geometry might give evidence 
for the localization of stimulated emission, before spon- 
taneous emission sets in. 

To prove the duality relation ([?]) we consider the prop- 
agation of monochromatic radiation (scalar amplitude E, 
wavenumber k), described by the Hclmholtz equation 



HE(r) = 0, 



n 



k 2 e{f). 



(4) 



(We suppose that all polarization-sensitive phenomena 
are absent.) Disorder leads to spatial fluctuations of the 
real part e'(r) of the dielectric constant. In the absence 
of disorder e' = 1. A non-zero imaginary part e" makes 
the system non-conservative. For simplicity we assume 
a homogeneous e". Its sign determines whether the sys- 
tem is absorbing (e" > 0) or amplifying (e" < 0). The 
parameter a introduced above is related to s" by 



a = -2fclm VT+ie 7, 



(5) 



where the argument of the square root is chosen in the 
interval (-tt/2, 7t/2). For \e"\ < 1 one has a = -ke". 

The dual symmetry underlying Eq. (|^) is formulated 
in its general form in terms of scattering matrices. We 
assume that the system consists of a scattering region 
of length L, in which e = e'{r) + is" , embedded in an 
iV-mode waveguide with e = 1 (see Fig. [I], inset). The 
scattering matrix S is a 2N x 27V matrix relating incom- 
ing and outgoing modes. It has the block structure 



S = 



r t' 
t r' 



(6) 



where r, r' are the reflection matrices and t, t' the trans- 
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FIG. 1. Numerical simulation of the transmittance of a 
single-mode waveguide (W — d, k 2 = 0.5 d~ 2 , Se = 0.2, 
I = 521 d, N = 1), averaged over 10 4 realizations of the disor- 
der. The right half of the figure is for amplification (circles: 
7 = 0.1; triangles: 7 = 0.2), the left half is for absorption 
(crosses: 7 = —0.1; squares: 7 = —0.2). The solid lines 
are the analytical asymptotes from Eq. ([15|). Their slope is 
independent of the sign of 7, in agreement with the duality 
relation (0). The inset shows the geometry considered. 



mission matrices, 
arc defined as 



The transmittances and reflectances 



T = iV" 1 Trtt t , 
T' = AT 1 THY* 



R = N- lr Tirr\ (7a) 
R' = AT^Tr/A (7b) 



Here T and R are the transmitted and reflected flux di- 
vided by the incident flux from the left. Similarly, X" and 
R' correspond to incident flux from the right. By taking 
the trace in Eq. (Q) we are assuming diffuse illumination, 
i.e. that the incident flux is equally distributed over the 
N modes. In the absence of gain or loss (cr = 0) the 
scattering matrix is unitary, SS^ = 1. This relation ex- 
presses flux conservation and relies upon Hermiticity of 
the Helmholtz operator, TL = TV at a = 0. For non-zero 
a we have Tt(a) = 1-0 {— a), which implies the duality 
relation 



S(a)S\-a) = 1. 



(8) 



Let us now examine the consequences of the duality 
relation (|J) for the reflection and transmission matrices 
of two systems which differ only in the sign of a. (We 
call these systems "dual".) Take a > 0, so that —a cor- 
responds to absorption. In the limit L — > 00, all elements 
of t— a and t'_ a tend to zero, while r_ CT and r'_ a remain fi- 
nite. Expanding the inverse of S(— cr) to first order in the 
transmission matrices and equating the result to S'{cr), 
we find 



tt 



0{t 2 ), 



0{t 2 



1 -a L -o' -cr 



(D{t 2 ). 



(9) 
(10) 



We introduce the transmission and reflection eigenvalues 
T n , 7^, 7l n , 1Z' n , being the eigenvalues of, respectively, 
T = tit, t' = t'i't, R = rr\ R' = r'r't. Because of 
time-reversal symmetry S(cr)S*{— cr) = 1. Together with 
Eq. (j^) this implies that S is a symmetric matrix. It 
follows that t' = t T , hence T n = T„' and T = T . The 
reflectances R and R! may differ. Eq. (^) directly yields a 
duality relation for the reflection eigenvalues in the limit 
L — > 00, 



(11) 



Eqs. @ and jlfj| ) together imply that the matrices 
TLjRI* and T^R^T 1 have the same eigenvalues. The 
duality relation for the transmission eigenvalues follows 
from the following lemma: 

Let A(L) be a matrix function of L with exponentially 
decreasing eigenvalues a n (L). The eigenvalue localization 
lengths £ n are defined by £-* = — liniL^oo In a n (L). 
Let B(L) be another non-singular matrix function whose 
elements remain finite as L — > 00. Then the matrix AB 
has the same eigenvalue localization lengths as A. 

It follows that the matrices T a , T^R^T 1 , T'_ ff Rl^, 
T'_ a> and hence T_ ff all have the same eigenvalue local- 
ization lengths. Explicitly, 



lira L 1 In T n (cr) 



lim L- l \a.T n {-a). (12) 



The transmittance T = N^ 1 J2 n 1S dominated by the 
largest transmission eigenvalue, which is the T n with the 
largest localization length: £ = max(£i, £2, • ■ ■ , £,n)- This 
completes the proof of Eq. (||) , since we have shown that 
all, and in particular the largest, transmission eigenvalues 
of dual systems have the same localization length. 

The case N = 1 of a single-mode waveguide can be ana- 
lyzed in more detail. We assume that the wavelength A is 
much smaller than both I and l/|cr|. (This is not a restric- 
tive assumption for an optical system.) The joint proba- 
bility distribution P(R, T, L) of reflectance and transmit- 
tance evolves with increasing L according to a Fokker- 
Planck equation, 



dP 



d_ 

dR 



[{1-R) 2 + 2 1 R]P- 



d 2 
dR 2 



R(l - R) 2 P 



-^(7-1 



i) 2 



R)P+^T 2 RP 



if 



dRdT 



TR(1 - R)P, 



(13) 



where we have abbreviated 7 = al. For 7 < (absorp- 
tion), this equation is equivalent to the moment equation 
of Freilikher, Pustilnik, and Yurkevich R], For 7 > 
(amplification) their method of moments cannot be used, 
because all moments of R diverge as L — > 00. The deriva- 
tion of Eq. ( p"3| ) proceeds along the lines of Ref. [ p^[ , 
where the case 7 = was considered. On integration over 



2 



T it reduces to a well known |14-|l6| Fokker-Planck equa- 
tion for P(R,L) = j dT P(R, T, L). The limit L -> oo 
of P(R,L) was studied in Refs. [p~5|— p~T[ . In terms of the 
variable fi = \/(R— 1) it reads 



2 7 e- 2w 6»(Ai) 
_2 7e -27(i+M) 



-1 - M ) 



for 7 > 0, 
for 7 < 0, 



(14) 



where 6(x) = 1 for x > and otherwise. Using this 
asymptotic distribution we have computed from Eq. fll3| ) 
the first two moments of InT in the large-L limit. The 
result for the average is 



(\nT)=-(l + \ 7 \)L/l + 2c( 7 ) 




(15a) 



c(7) _ ( C + In 2 7 - e 2 ^Ei(-2 7 ) for 7 > o] (15b) 
where C is Euler's constant and Ei(x) = dt e* jt is 



the exponential integral (0(7) s» — 2 7 ln 7 if < 7 <C 1). 
For 7 < 0, Eq. @ reduces to the result of Refs. §,@. 
The result for 7 > is new, and demonstrates that the 
inverse localization length £ _1 ((t) = (1 + I7I) l^ 1 = I' 1 + 
\<j\ is indeed independent of the sign of a — in accordance 
with the general duality relation (||). The result for the 
variance is 



varlnT = 2 



l + 2|7|e 2 ^Ei(-2| 7 |)U/Z + 0(l), (16) 



in agreement with Ref. jjj for 7 < 0. Note that 
\/var In T < (InT) for L/l > 1. Evaluation of higher 
moments shows that the distribution of In T tends to a 
Gaussian for L — > 00. (The tails are non-Gaussian, but 
contain negligible weight.) 

These results hold in the large-L limit. For short 
waveguides instead of Eq. ( plf ) one has (InT) = — (1 — 
7)L/Z. The crossover length can be estimated as L c ~ 
ic(7)/|7|. Below this length stimulated emission en- 
hances transmission through the waveguide. On larger 
length scales stimulated emission reduces transmission. 
In contrast, the reflectance is enhanced on every length 
scale IJHL7). 

To test these analytical predictions for N = 1, and 
to investigate also the multi-mode case, we have numeri- 
cally solved a discretized version of the Helmholtz equa- 
tion (^J) , on a two-dimensional square lattice (lattice con- 
stant d, length L, width W). The real part e' of the di- 
electric constant was chosen randomly from site to site 
with a uniform distribution between 1 ± 5s. The scatter- 
ing matrix for the multi-mode case was computed using 
the recursive Green's function technique, originally de- 
veloped for the electronic Anderson model [Q . (For the 
case N = 1 a transfer-matrix method || turned out more 
convenient.) Simulations with e" = were used to obtain 
I, from the relation |lj| 



The parameter a was determined from Eq. (|). Results 
for the single-mode case are shown in Figs. |l and ^, and 
for the multi-mode case in Fig. [|. The duality relation 
between the localization lengths for absorption and am- 
plification is verified with good accuracy, both for the 
single- and for the multi-mode case. Furthermore, for 
N = 1 we find good agreement with the results (15)-(|l6|) 
of the Fokker-Planck equation. 
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FIG. 2. Probability distribution of the logarithm of the 
transmittance of a single-mode waveguide, for L/l — 15.4 
and 7 = 0.2 (triangles, dashed curve), 7 = —0.2 (squares, 
solid curve). The data points are provided by a numerical 
simulation (same parameters as in Fig. |l|), the curves are a 
Gaussian distribution of In T with mean and variance given by 
Eqs. @ and ©. There is a slight offset between the distri- 
butions for absorption and amplification because the system 
is not fully in the large-L limit. 

In conclusion, we have demonstrated that stimulated 
emission of radiation in a disordered waveguide reduces 
the decay length, in the same way as absorption does. 
This paradoxical result is an immediate consequence of 
the exact duality relation (^) between the scattering ma- 
trices of two systems with complex conjugated dielectric 
constants. The dual symmetry between absorption and 
amplification has been supported by an explicit compu- 
tation of the decay lengths, both analytically (for the 
single-mode case) and numerically (for the single- and 
multi-mode cases). 

We acknowledge useful discussions with P. W. Brouwer 
and K. M. Frahm. This work was supported by the Dutch 
Science Foundation NWO/FOM. 

Note added: Numerical results for the single-mode case 
supporting Eq. (|j) have been published by Z.Q. Zhang, 
Phys. Rev. B 52, 7960 (1995). 



lim L _1 (lnT) = [h(N + l)l] 

L— ^ao 



(17) 



3 







-5 



o 

■U/ 

+ 

-5-10 



15 



absorption 



□ 



□ 



f □ 2.0 
□ 

□ 1.5 
^1.0 
0.5 
0.0 



amplification 



A O r 



5 a a 



20 10 10 20 
^L/l L/l^ 

FIG. 3. Numerical simulation of the transmittance of a 
multi-mode waveguide (W = 25 d, k 2 = 2.0 d~ 2 , Se = 0.375, 
I — 29.6 d, N = 12), averaged over 50 realizations of the dis- 
order. The parameter £o = {N + 1)1/2 is the localization 
length of the system in the absence of absorption or amplifi- 
cation. The right half of the figure is for amplification (cir- 
cles: a = 0.0035 d" 1 ; triangles: a = 0.0071 d" 1 ), the left 
half is for absorption (crosses: a = —0.0035 d~ ; squares: 
a = —0.0071 d -1 ). The inset shows the eigenvalue localiza- 
tion lengths, = — lirriL^oo L" 1 In 7^. These lengths £ n 
were computed from the L-dependence of T n for L up to 40 1 
and a single realization of the disorder. The duality between 
absorption and amplification is verified with good accuracy. 
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